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Abstract 

(N ■ 

>■ ' The OPE treatment that has been so successful in describing inclusive B — > lvX c decays 

yields sum rules (in particular the Uraltsev sum rule and its higher moments) implying 
the dominance of the P wave j q = 3/2 charm states in X c over their j q = 1/2 counter- 
parts. This prediction is supported by other general arguments as well as quark model 
calculations, which illustrate the OPE results, and by preliminary lattice findings. Its fail- 
ure would indicate a significant limitation in our theoretical understanding of B — >• lvX c . 
Some experimental issues have been clarified since a preliminary version of this note had 
appeared, yet the verdict on the composition of the final states beyond D, D* and the two 
narrow j q = 3/2 resonances remains unsettled. Establishing which hadronic configura- 
tions - D/D* + 7r, D/D* + 2tt, ... - contribute, what their quantum numbers are and their 
mass distributions will require considerable experimental effort. We explain the theoreti- 
cal issues involved and why a better understanding of them will be of considerable value. 
Having significant contributions from a mass continuum distribution below 2.5 GeV raises 
serious theoretical questions for which we have no good answer. Two lists are given, one 
with measurements that need to be done and one with items of theoretical homework. 
Some of the latter can be done by employing existing theoretical tools, whereas others 
need new ideas. 
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1 Outline 

Both our theoretical and experimental knowledge on semileptonic B decays have advanced 
considerably over the last 15 years. This progress can be illustrated most strikingly by 
the recent success in extracting the value of |V(c6)| with better than 2 % accuracy from 
measurements of inclusive B — > lvX c transitions [2, 3, 4]. At the same time some potential 



problem of a rather subtle nature have emerged. One concerns BK(B — > lvX c ), which 
has been very well measured now [5]: 

BR(B d — > l~9X c ) = (10.33 ±0.28)% (1) 
BR(B U -> rvX c ) = (10.99 ±0.28)% (2) 



While the ratio of these branching fractions is well understood in terms of the B d — B~ 
lifetime difference, their absolute scale falls below early predictions inferred from the 
Heavy Quark Expansion (HQE). Yet those were based on values for the charm quark 
mass that appear too heavy now. Using smaller values of m c (together with a more 
careful definition of heavy quark masses) and including some novel radiative corrections 
[7] enhances the rate for b — > ccs and thus lowers BR SL (B). It also enhances the charm 
content in the final state of B decays - in agreement with the data [8]. A conclusive 
theoretical analysis of BRsl(B) has not been performed yet although the tools exist. Yet 
we do not suspect this observable to represent a real problem for theory. 

In this note we want to focus on the composition of the hadronic final state in semilep- 
tonic B decays beyond B — > WD/D*. Understanding the nature of the hadronic system in 
the final state - its quantum numbers as well as mass distributions - is important, since 
well grounded theoretical expectations and predictions can and have been given on these 
issues. Heavy quark symmetry that becomes an exact symmetry of QCD in the limit 
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vtlq — » oo provides at least a convenient classification scheme. The S wave configurations 
D and D* represent the ground states to be followed by four P wave [cq] excitations. In 
two of those the light degrees of freedom carry angular momentum j q = 3/2 resulting in 
narrow resonances with spin 1 and 2. For the other two one has j q = 1/2 leading to broad 
resonances with spin and 1. Our theoretical understanding of semileptonic B decays 
tells us rather unequivocally that the j q = 3/2 states should be more abundant in the 
final states than their j q = 1/2 counterparts. This prediction appears to be at variance 
with some data. We refer to this apparent conflict as the '1/2 vs. 3/2 puzzle' [1]. One 
also expects smallish contributions from other hadronic configurations. There still seems 
to be some tension between data sets concerning the mass distribution of such additional 
contributions and their decay patterns. The aim of this note is to explain in a concise way 
the arguments involved in deriving the theoretical expectations and the consequences of 
their possible failure to stimulate further experimental as well as theoretical studies. The 
issues had caught our attention and led to the original 'Memorino' (= short memo) more 
than a year ago [9]. While more data have been obtained since, and we have pondered the 
issues further, we find the problems now even more intriguing and in need of a resolution. 
The latter has to be driven by even more detailed analyses. A more appropriate name for 
this paper might now be 'Memorone' (= long memo) - alas we decided to stick with the 
original moniker. 

After giving an overview of the experimental situation for B — > luX c in Sect. 2 we 
marshall the theoretical arsenal for treating those decays: the operator product expansion 
(OPE) in Sect. 3, the BT model in Sect. 4 and lattice QCD in Sect. 5 before adding other 
general arguments in Sect. 6; in Sect. 7 we undertake a more detailed comparison of the 
theoretical predictions and expectations with the existing data on5^ lVD^ + 7r's from 
ALEPH, BaBar, BELLE, CDF, DELPHI and DO; in Sect. 8 we comment on corresponding 
expectations for nonleptonic B decays; finally in Sect. 9 we list needed homework for both 
theorists and experimentalists. We aim at being as concise as reasonably possible, while 
providing a guide through the literature for the more committed reader. 



2 The data 

About three quarters of the inclusive semileptonic B width are made up by the two 
channels B — > IvD/D*, for which a recent BaBar analysis finds [14]: 

r(B~ -> l~DD) T(B- -> l~vD*) 

0.227 ± 0.014 ± 0.016, — y — = 0.582 ± 0.018 ± 0.030 



T(B~ -> l-VDX) ' ' ' ' T(B~ -> 1-vDX) 

T(B d -> l~PD) T(B d -> l~vD*) 

' = 0.215 ± 0.016 ± 0.013, — ^ — = 0.537 ±0.031 ±0.036 



V{B d ^l-DDX) ' Y{B d — > l~VDX) 

(One expects T(B — > l~VDX) to saturate T(B d — > l~i?X c ) for all practical purposes, and 
this is completely consistent with observation.) This large dominance of D and D* final 
states, the ground states of heavy quark symmetry, represents actually the most direct 
evidence that charm quarks act basically like heavy quarks in B decays. This can be 
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invoked to justify using the heavy quark classification already for charm and applying 
arguments based on the SV limit [15]. 



For the remainder we have 

T(B- -> l~VD**) 



T(B- -> 1-uDX) 
T(B d -> I'vD**) 
T(B d -> 1-vDX) 



= 0.191 ±0.013 ±0.019 (3) 
= 0.248 ± 0.032 ± 0.030 , (4) 



where D** denotes any D^nn combination that is not a D*. The apparent sizable dif- 
ference in the central values of these two ratios is in conflict with theoretical expectations 
based on isospin symmetry. For the latter imposes practically equal rates for the cor- 
responding semileptonic channels of B~ and B d . We view this difference as due to a 
statistical fluctuation or a systematic bias. Averaging over B~ and B d rates yields [14]: 

T(B -> I'PD) F(B -> l~DD*) 

- 0.221 ±0.012 ±0.006 , — L —4 = 0.572 ± 0.017 ± 0.016 



T(B -> 1-uDX) ' ' ' ' T(B -> 1-uDX) 

Y(B -> Z-PD**) 



r(s -> i-vDX) 



0.197 ±0.013 ±0.013 (5) 



BELLE finds completely consistent branching ratios for 5 _ — > l~v>D° / D*° and -B^ — > 
l~j)D + . Yet their number for — > l~j)D* + appears on the low side compared with the 
expectation based on isospin. As before we take that as a sign that the data are not 
completely mature yet rather than as a real effect. 

The numbers in Eq.(5) are consistent with the HFAG averages [5]: 

BK{B d -> l~VD + ) = (2.08 ± 0.18)% , BK(B d — ► l~DD* + ) = (5.29 ± 0.19)% 

BR(_B d — > l~i>D** + ) = (2.8 ±0.3)%, (6) 

where the D** rate is inferred by subtracting the Z) and D* rates from the total semilep- 
tonic rate. Another even more recent BABAR analysis [6] finds numbers manifestly 
consistent with isospin invariance: 

(7) 
(8) 
(9) 
(10) 

The ratio of corresponding branching ratios is fully consistent with the B~-B d lifetime 
ratio. 

3/2 1/2 

The four P wave excitations D 12 and D 01 are obvious candidates to provide D** 
contributions. ALEPH [28] has reconstructed D** states decaying into D^*- ) 7r ± . They did 
not observe a significant excess of events over the expected background in D^ + n + or D°ir~ 
combinations (called 'wrong sign'). From the measured rate of 'right sign' combinations 



BR(B- 


-> 1-PD°) 


= (2.33 ± 0.09 stat ± 0.09 syst ) 


% 


BR(B- - 


-> rpD*°) 


= (5.83 ±0.15^^0.30^) 


% 


BR(B d - 


-> rvD + ) 


= (2.21 ± 0.1 ±0.12 syst ) 


% 


BR(B d - 


> l~uD* + ) 


= (5.49 ± 0.16 stat ± 0.25 syst ) 


% 
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and assuming that only D** decaying to D^n contribute (to correct for channels with a 
missing vr°) they get (with Prob(6 -> B) = (39.7 ± 1.0)%): 

BR(B -> WD** -> WD^n) = (2.2 ± 0.3 ± 0.3)% (11) 

Assuming the D^ 2 state to decay only into D*n they find among their D** sample: 

BR(B ^WD\ /2 ) = (0.70 ±0.15)% (12) 
BR(B^WD 3 2 /2 ) < 0.2%. (13) 

DELPHI has published a re-analysis of their data [29] superseding their previous study 
[30]. Assuming only D^n to contribute (to correct for channels with a missing ir°) they 
obtain 

BR(B -> WD** -> WD^n) = (2.7 ± 0.7 ± 0.2)% (14) 
with clear evidence for two narrow states tentatively identified with D 3 ^ 2 : 

BR(B -> /z/Di 72 ) = (0.56 ± 0.10)% , BR(B -> ZP£>^ /2 ) = (0.30 ± 0.08) , (15) 

The DO collaboration has measured production rates of narrow D** states in the decay 
B -> pi-9D*7r. Assuming BR(L>i /2 -> £>*7r) = 100% and BR(D 3 2 /2 -> £>*tt) = (30 ± 6)% 
they obtain [31] 

BR(5 -> /i-zaDi /2 ) = (0.33 ± 0.06)% , BR(B -> n~vD\J 2 ) = (0.44 ± 0.16)% . (16) 
These three sets of data agree in pointing to 

BR(5 -> rvD\%) ~ 0.8 - 1% , (17) 

yet do not paint a clear picture on the relative strength of B — > \~vD\J 2 vs. 5 — > l~vD^ 2 . 

The observed rates for 5 -> WuD/D*/D^ 2 are in pleasing agreement with the the- 
oretical predictions described below. The fact that they do not quite saturate Fsl(B) is 
not surprising. Even so one would like to learn from the data what these additional final 
states are, in particular what their distribution in mass is. If they populate the range 
above 2.5 GeV, then we have at least candidates for them. If however they fall mostly 
below 2.5 GeV, then we can come up at best with rather exotic explanations. These 
are subtle questions concerning smallish rates. Yet answering them will teach us impor- 
tant and potentially very surprising lessons on nonperturbative dynamics. This will be 
explained in Sect. 7 after having marshaled the theoretical predictions and expectations. 



3 The OPE treatment 

While the OPE allows to describe inclusive transitions, no systematic extension to ex- 
clusive modes has been given so far. Yet even so, the OPE allows to place important 
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constraints on some exclusive rates: B — > lisD*, IvD (the latter involving the 'BPS' ap- 
proximation on which we comment later) are the most topical and elaborated examples 
[10,11]. 

OPE results can be given also for subclasses of inclusive transitions due to various sum 
rules [12, 13] that can genuinely be derived from QCD; hence one can infer constraints on 
certain exclusive contributions. Those can be formulated most concisely when one adopts 
the heavy quark symmetry classification scheme also for the charm system in the final 
state of semileptonic B decays. 

In the limit uiq — > oo one has heavy quark symmetry controlling the spectroscopy for 
mesons as follows: The heavy quark spin decouples from the dynamics, and the hadrons 
can be labeled by their total spin S together with the angular momentum j q carried by 
the light degrees of freedom, namely the light quarks and the gluons. The pseudoscalar 
and vector mesons D and D* then form the ground states of heavy quark symmetry in 
the charm sector with j q = 1/2. The first excited states are four P wave configurations, 
namely two with j q = 3/2 and S — 2, 1 - D% 2 , D^ 2 - and two with j q = 1/2 and 
S — 1,0 - D^ 2 , D l J 2 ] the two 3/2 states are narrow resonances and the two 1/2 states 
wide ones. Then there are higher states still, namely radial excitations and higher orbital 
states; furthermore there are charm final states that cannot be properly called a hadronic 
resonance, but are D/D* combinations with any number of pions etc. carrying any allowed 
J PC quantum numbers. 

The P wave states D^ 2 and D\q are obvious candidates for D**, but they should not 
saturate it completely. One actually expects QCD radiative corrections to populate the 
higher hadronic mass region above the prominent resonances through a smooth spectrum 
dual to a superposition of broad resonances. 

Up to small isospin breaking effects one predicts these semileptonic rates to be the 
same. The usual Isgur-Wise function £(u>) is the core element in describing B — > IvD/D*. 
It can be generalized to describe also the production of excited charm final states in 
semileptonic B decays: r i/2[3/2] ( u '™) with w n = vb • v D ( n ) is the amplitude for B — > 



where -D 1 /2[3/2] denotes a hadronic system with open charm carrying j q 



1/2 [3/2] and label n\ it does not need to be a bona fide resonance. 

Various sum rules can be derived from QCD proper relating the moduli of these 
amplitudes and powers of the excitation energies = M D („) — Mc to heavy quark 
parameters. Adopting the so-called 'kinetic scheme', as we will throughout this note, one 
obtains in particular [12]: 



= -E 



+ EM/ m 2 ) (i) 



e< m )<^t 



= E (4/ 2 ) |t$(1)| +2 E (4?2 } ) |r$>(l) 

n m 

= -2 e ($) 2 h$(i)f+2 e" ^TMr^ 1 ) 



(18) 
(19) 
(20) 
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where the summations go over all hadronic systems with excitation energies e( n,m ) < /i 
1 . The sum rules relate the observables a4(a*) and /^(aO that describe the fully inclusive 
B — > l~i?X c transitions with sums over subclasses of exclusive modes. For the purpose of 
these sum rules it is irrelevant whether these hadronic systems are bona fide resonances 
or not; what matters is that they are P wave configurations with j q = 3/2 or 1/2. 

These sum rules allow us to make both general qualitative as well as (semi) quantitative 
statements. On the qualitative level we learn unequivocally that the '3/2' transitions have 
to dominate over the '1/2' ones, as can be read off from Eqs.(18,20). Furthermore we know 
that 

n 2 M > nUv) (2i) 

has to hold for any /i [15], as is obvious from Eqs.(19,20). On the quantitative level it is 
not a priori clear, at which scale \i these sum rules are saturated and by which kind of 
states. We will address these issues below. 

We have learnt a lot about the numerical values of the heavy quark parameters: the 
most accurate value for the chromomagnetic moment n G can be deduced from the B* — B 
hyperfine mass splitting: 

/4(1 GeV) = (0.35 ± 0.03) GeV 2 (22) 

The analyses of Refs.[2, 3, 4] based on comprehensive study of energy and hadronic 
mass moments in B — > luX c yield the following values: 

0.297 ± 0.024| ea;p ± 0.046|hq£ GeV 2 [2] 
/4(lGeV) = { 0.358 ± 0.060| /it ± 0.003| fos GeV 2 [3] (23) 
0.330 ± 0.042| ea;p ± 0.043U eo GeV 2 [4] 

0.401 ± 0.019U P ± 0.035|^Q£ GeV 2 [2] 
/i 2 (lGeV) = { 0.557 ± 0.091| /it ± 0.013| fas GeV 2 [3] (24) 
0.471 ± 0.034| e:rp ± 0.062| t?ieo GeV 2 [4] 

These experimental numbers are consistent with each other and with the predictions 
of Eqs. (22,21) without the latter having been imposed. All three determinations of 
/x 2 (l GeV) are within one sigma of 0.45 GeV 2 , which we will use as a reference point 
for our subsequent considerations: 

/i 2 (l GeV) | rcf = 0.45 GeV 2 . (25) 

For the lowest excitation energies we have 

~ 450 MeV , ej ^ ~ (300 - 500) MeV ; (26) 

these values for efj 2 allow for '1/2' states to be both lighter and heavier than the narrow 
'3/2' states. 



1 Thc sum rule of Eq.(18) does not require a cut-off or normalization scale /i, as is already implied by 
its left hand side [12]. 
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From Eqs. (19,20) one obtains 

a4(a*) 



6 E ( 



n 

(™)</i 



e (n) V 

3/2 J 



_(») 
r 3/2 



(l)f-^(^)-^)) (27) 



^)-^) = 9 £ (4/0 |^(1) 



(28) 



as convenient expressions to read off natural 'scenarios' for the implementation of the 
OPE description and its sum rules. 

One can reasonably assume these sum rules to be saturated approximately by the 
lowest states n = for /i < 1 GeV. This rule of thumb (not to be confused with sum 
rules) is based on general experience with sum rules and on considerations of how 
and //^(/-O vary with the scale [i. It does not mean that the various sum rules would 
saturate to the same degree at a given for they reflect different dynamical situations. 

Using Eq.(25) we then infer from Eq.(27) 



0.6 



(29) 



and from Eq.(28) 

7-^(1) < 0.14 - 0.32 , (30) 

which are reasonable numbers as our subsequent considerations will illustrate. For /x 2 (1 GeV) = 
0.4 GeV 2 one has t[°\(1) < 0.1 - 0.2 and r 3 ( ^(l) ~ 0.56. On the other hand /4(1 GeV) = 

0.55 allows for a sizeable production rate of the lowest '1/2' state - Ty\(l) < 0.2 — 0.45 - 



although does not enforce it, and r. 



(0) 
3/2 V 



0.63. These numbers for t^\(1) should also be 



seen more like an upper bound, since the assumed saturation of the sum rules can hardly 
be exact. 

Imposing also Uraltsev's sum rule, Eq.(18), with (approximate) saturation assumed 



for the lowest states 



T- 



(0), 
1/2 1 



1) 



t$(1)| S r~0.6, t$(1)|sr~0.32; 



r (o) 

r 3/2 [ 



1) 



— | J leads to 

(0), 



(31) 



i.e., very much the upper end of Eq. (30) inferred from /i 2 (l GeV) = 0.45 GeV 2 and 
pointing to efj 2 ~ 300 MeV, i.e. a relatively low mass for the '1/2' states. A scenario 
with a low /i 2 (l GeV) =0.4 GeV 2 is however hardly compatible with it, whereas a high 
GeV) = 0.55 GeV 2 can be accommodated with ef] 2 ~ 450 MeV. 

While these numbers are reasonable, one cannot rule out significant contributions from 
higher states like the first radial excitations of the P wave states. We will address such a 
scenario in Sect. 4. 

In summary: 

• The OPE treatment leads to the general prediction that among P wave configura- 
tions production of '3/2' states dominates over that of '1/2' states in semileptonic 
B decays and that the former yields a significant contribution. 
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• Approximate saturation of the sum rules by the P wave states represents a scenario 
consistent with the data and leads to semi- quantitative estimates for the degree of 
the dominance of '3/2' over '1/2' production. 

• At the same time there is no reason from the OPE to expect that D, D* and 
the two narrow D 3 / 2 states saturate the semileptonic width. We know that the 
OPE treatment describes the hadronic mass moments in semileptonic B decays 
very sucessfully [2, 3, 4], and the data clearly show BR(_B — > l~vD x ) ~ 1 — 2% 
with D x 7^ (D, D*, D^^). Already on general grounds one would not expect 
Dx to be mostly a narrow resonance. More specifically one can start from the 
numbers in Eq. (6) and split the B — > l~vD** contribution into two components 
with BR(B -> l-DDWV) = 0.8% and BR(B -> l~vD x ) = 2%. From the hadronic 
mass moments determined in the kinetic scheme [2, 3] one can then infer by matching 
what the hadronic mass moments for the Dx contribution have to be. One typically 
finds (M(D X )) ~ 2.4 - 2.6 GeV with a spread of about 200 MeV. DELPHI has 
inferred the mass moments of the D** contributions and found a central value of 2.5 
GeV with a spread of 230 MeV. 

Then the question arises what makes up this Dx contribution. Some broad '3/2' 
configuration presumably of a non-renonant nature? Or states that do not con- 
tribute to the sum rules like J p = 0~, 1~ states? Radial excitations would fit this 
bill, yet run counter to arguments to be discussed in Sect. 7. 

• No reliable prediction on their decay patterns - i.e. whether they yield D^ir or 
D^^tt etc. - can be inferred from the OPE treatment per se. 



Based on the OPE treatment alone one cannot be more specific numerically. To go further 
one relies on quark models for guidance. The dominance of the '3/2' over the '1/2' states 
emerges naturally in all quark models obeying known constraints from QCD as well as 
Lorentz covariance. This can be demonstrated explicitly with the Bakamjian-Thomas 
covariant quark model [17], which satisfies heavy quark symmetry and the Bjorken as 
well as spin sum rules referred to above. It allows to determine the masses of various 
charm excitations and to compute the production rates in semileptonic [18, 19, 20] as well 
as nonleptonic B decays [19]. The BT model provides a quantitative illustration of the 
heavy quark limit, in particular concerning the sum rule of Eq.(18). One finds 



together with predictions for the slopes. These values are fully consistent with the esti- 
mates given above. For the semileptonic modes the BT model yields: 



4 The BT model 




(32) 
(33) 



BR(B -> WD) 



(1.95 ±0.10)% , 



(34) 
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BR(B -> WD*) = (5.90 ±0.20)% (35) 

BR(B^WD 3 2 /2 ) = (0.631°;*)% (36) 

BR(B^WDl /2 ) = (0.401°;^)% (37) 

BR(B ^lvD\' 2 ) = (0.06 ±0.02)% (38) 

BR(B -> WDl /2 ) = (0.06 ±0.02)% (39) 



The following basis and features of these predictions should be noted: 

• The predictions for B — > IvD/D* are based on the following parametrization of the 
Isgur-Wise function: £(w) = (^y) \ where p 2 denotes its slope. For the latter we 
have used the value from Ref. [23]. The predictions agree with the data. 

• The branching ratios for the P wave states and the theoretical uncertainties are 
obtained by using the BT model value for the form factors for w — 1 and allowing 
for a ±50 % variation in the slope given by the BT model. This is the origin of the 
large relative errors in the predicted rates. We find a strong dominance of '3/2' over 
'1/2' production for the P wave states as inferred already from the Sum Rules. 

• These branching ratios add up to 9.00 ± 0.40 % and thus fall short of saturating 
the observed Fsl(B), Eq. (2). Such a deficit is not surprising as argued before. 
The more specific question is whether the BT model can account for the required 
additional width. The answer is not known yet. What can be said is that the 
hadronic mass for these extra contributions cannot be much lower than 2.6 GeV 
and that a sizable number of channels might be involved. 

• The model as it is does not allow to compute I/tuq corrections; i.e., effectively it 
treats the rriQ — > oo limit, since only then it is covariant. 

With an explicit quark model one can address also higher states. The BT model finds 
for the first radial excitations of the '3/2' and '1/2' P wave states still sizable amplitudes 
T3J2 — 0.21 and Ty 2 ~ 0.20. While these states enhance GeV) significantly, their 
contributions to Tsl(B) are still found to be insignificant. 

5 Lattice QCD 

In principle the two form factors Ti/ 2 (1) and r 3 / 2 (l) can be computed in a straightforward 
way using the HQET equation of motion (v • D) h v — [24]: 

v ^ + \h vl l r>D^h v \0-) v = i^T 1/2 (l)(A 0+ -A -), 
v ( 2 +\hjll p_Z \h v \0-) v = -zv / 3e^r 3/2 (l)(A 2+ -A -), (40) 
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where v = (1, 0) is the heavy quark velocity, e* the polarization tensor of the 2 + state and 
Ajp the dominant term in the OPE expression for the J p heavy-light meson binding en- 
ergy. On the lattice the covariant derivative Di applied to the static quark field h(x, t) is 
expressed as Dih(x,t) — > ^ (Ui(x,t)h(x + i,t) — U}(x — i,t)h(x — i,tfj; Ui(x,t) denotes 
the gauge link. One calculates as usual the two-point functions Cj P (t) = (0|Ojp(£)Ot P (0)|0), 
the three-point functions C^ -(ti,t 2 ) = (0|Ojp(t 2 )O r (ti)Oj_(0)|0) and (J p \O r \0~) oc 

Alas, numerical complications appear, because orbital as well as radial excitations can 
contribute. To extract properly the matrix element for the P wave state (J p |O r |(T), 
one has to effectively suppress the coupling of radial excitations (with quantum numbers 
n > 1, J p ) to the vacuum. This can be achieved by choosing an appropriate interpolating 
field Ojp such that (n > 1 J p |Ojp|0) = holds or by having huge statistics to diminish 
statistical fluctuations at large times (where the fundamental state is no more contam- 
inated by radial excitations). This poses a problem in particular for the 2 + state, for 
which the usual interpolating field seems to couple also the first radial excitation quite 
strongly to the vacuum. Moreover reaching the required stability of R(ti,t 2 ) as a function 
of t 2 poses a serious challenge. Hopefully all to all propagators technology will be of great 
help, as it has already proved to be in studies of the static-light spectrum [25] and in the 
determination of hadronic matrix elements [26] . 

We will need very careful and dedicated lattice studies to obtain meaningful and 
reliable results for r 3 / 2 ,i/2- As an already highly relevant intermediate step one can con- 
centrate first on Ti/2 to see whether lattice QCD confirms its suppression as inferred from 
both the sum rules and the BT model. A preliminary study in the quenched approxima- 
tion with j3 = 6.0 (a^ 1 = 2 GeV -1 ) and m q ~ m s yields [27]: 

t$(1)|lqcd ~ 0.41 ±0.05 (41) 

3/2V 



^3 ( /2(1)Ilqcd ~ 0.57 ±0.10, (42) 



where only the statistical errors are given. Again we find dominance of the '3/2' over the 
'1/2' amplitude even in numerical agreement with the values inferred from the sum rules, 
see Eq.(31); t^^IOIlqcd appears significantly larger than the BT estimate in Eq.(33). 

Apart from unquenching and lowering the value of m q one can improve and refine this 
analysis also by simulating a non-static charm quark, i.e. applying HQET to the B meson 
only. This would allow to evaluate l/m c corrections. However the first improvement has 
to be done with considerable care: While unquenched simulations with light quark masses 
lower than m s /5 (i.e. a corresponding pseudoscalar Goldstone boson mass lower than 300 
MeV) have become customary now, one has to worry about a possible mixing between a 
D** resonance and a Dir state. 
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6 Two other general arguments on | ^ / 2 / ^"3/2 1 2 



The numerics of the theoretical predictions on semileptonic B decays given above have to 
be taken 'cum grano saris'. Yet their principal feature - the preponderance of '3/2' over 
'1/2' states - has to be taken very seriously, since they are a general consequence of the 
OPE treatment. It is further supported by two rather general observations that point in 
the same direction as the detailed theoretical considerations given above: 

• When interpreting data one should keep in mind that the contributions of D\p 
to T(B — > WD**) are suppressed relative to those from by a factor of two 
to three due to kinematics [18]. Thus one finds for reasonable values of Ty\ that 
F(B — > WD 1 / 2 ) falls below T(B — > WD 3 / 2 ) by one order of magnitude, as illustrated 
above, see Eqs. (38,39). For the two widths to become comparable, one would need 
a greatly enhanced Ty\. 

• There is a whole body of evidence showing that in so-called class I nonleptonic B 
decays like B^ — > D^+n - naive factorization provides a very decent description 
of the data. Invoking this ansatz also for Bd — > D** + n~ — > D^ix+n^ one infers 
from BELLE's data [21] that the production of '1/2' states appears to be strongly 
suppressed relative to that for '3/2' ones. It implies that IT1/2/T3/2I 2 is small and 
certainly less than unity. The same feature is found in more recent measurements 
from BABAR [22] 2 . This agrees with the theoretical expectations described before; 
more importantly it shows in a rather model independent way that there is no large 
unexpected enhancement of |ti/ 2 |. Those values also allow to saturate the sum rule 
of Eq.(18) within errors already with the n = states. 

The form factors are actually probed at w — 1.3 in this nonleptonic transition; yet a 
natural functional dependence on w supports this conclusion to hold for 1 < w < 1.3 
in semileptonic channels. 

These arguments are based on the heavy quark mass limit. The as yet unknown finite 
mass corrections could modify these conclusions somewhat. 



7 Detailed comparison with the data on semileptonic 

B decays 

Different experiments and theoretical treatments agree on 

• T(B — > WuX c ) being dominated by the two modes B — > WuD/D*, 

• with Z)( 3 / 2 ) final states providing about 10% to it and 

2 The comparison of the theoretical predictions with the measured rates is not straightforward, since 
the data are given as products of production and decays branching ratios, and one has to use the model 
to calculate both. The BABAR analysis also lumps together the decays of more than one state. 
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• BR(B — > / vX c ) ~ (1 — 2)% has to come from other hadronic configurations D X - 

The question still open concerns the nature of this last component D X - 

At first sight one might wonder why one should worry about the identification of 
channels that sum up to no more than 2% in overall branching ratio. Yet they constitute 
10 - 20% of all semileptonic transitions, and - maybe more importantly - theory makes 
quite non-trivial statements about them. We can learn important lessons about non- 
perturbative dynamics, even if those predictions are refuted by experiment. 

Theory makes a rather robust prediction that it cannot come from the broad P wave 
states D^^ 2 \ The OPE framework by itself can accommodate all three features listed 
above, as long as LK 1 / 2 ) is insignificant in the third item. It points to hadronic contri- 
butions that are broad in mass without a firm prediction on their average mass - both 
(M(D X )) < 2.4 GeV or > 2.5 GeV seem a priori feasible - or their decay patterns; i.e. 

£>(*) 

7T VS. 

£>(*) 

7T7T (VS. 

£>(*) 

7] etc.). 

Both ALEPH and DELPHI can account for all of T(B -> l~vX c ) with B -> l~DD/D* 
and B — ► l~vD** — > 1~uD^tt, see Eqs.(ll,14) with no established signal for D^nn 
states contributing. ALEPH places relatively tight bounds on higher combinations from 
the observed number of 'wrong sign' combinations: 

BR(B -> WD*mr) < 0.35% , BR(B -> WDirir) < 0.9% (90% C.L.) ; (43) 

DELPHI'S bounds are less tight: 

BR(B -> WD*7T7r) < 1.2% , BR(B -> ZzaDtttt) < 1.3% (44) 

Considering that 1 + can decay into Dim and analyzing the Dn mass distribution 
DELPHI fits a value of (19 ± 13)% for this component. In their analysis of hadronic mass 
moments such a possibility has been included with BR(_B — > IvDnir) = (0.36 ± 0.27)%. 
This turns out to be the dominant systematic uncertainty in their hadronic mass moment 
measurement. 

DELPHI found a significant rate for producing a broad hadronic mass distribution in 
BR(B -> ZPD-r) = (1.24 ± 0.25 ± 0.27)% , BR(B -> /z/D- ") = (0.65 ± 0.69)% . (45) 

1/2 

//the broad contributions were indeed to be identified with the D 10 as already implied 
in Eq.(45) - an a priori reasonable working hypothesis - one would have a clear cut and 
significant conflict with the OPE expectations as well as the numerically more specific BT 
model predictions, see Eqs.(38, 39). For DELPHI'S data would yield T(B — > luD 1 ^ 2 ) > 
F(B — > WD 3 / 2 ). This conflict has been referred to as the '1/2 > 3/2 puzzle' [1]. Since, as 
sketched before, the theoretical predictions are based on a rather solid foundation, they 
should not be discarded easily. Of course there is no proof that the broad D/D* + 7r 
systems are indeed the j q — 1/2 P wave states; they could be radial excitations or non- 
resonant combinations of undetermined quantum numbers. Thus the DELPHI data taken 
by themselves are not necessarily in conflict with theoretical expectations. 

However the plot thickens in several experimental as well as theoretical respects: 



13 



• In 2005 BELLE has presented an analysis of B — > WD/D*7i [32], which appears to 
be in conflict with previous findings. Reconstructing one B completely in T(4S') — > 
BB, they analyze the decays of the other beauty meson and obtain: 



BR(B- -> rPDir) = (0.81 ±0.18)% (46) 

BR(S~ — > 1~uD*-k) = (1.00 ±0.22)% (47) 

BR(B d -> I'DDtt) = (0.49 ±0.13)% (48) 

BR(B d — > 1~9D*-k) = (0.97 ±0.22)% (49) 



BELLE's separation of final states with D and D* is of significant value, since it 
provides an indirect and model dependent handle on '3/2' and '1/2' production. For 
with the help of a quark model one can calculate both the production rates for the 

(3/2) (1/2) 

D\ 2 and Dq { and their branching fractions into Dtt and D*ir. BELLE's numbers 
are actually quite consistent with the theoretical predictions the BT model yields 
for '3/2' P wave production. It is of course still desirable for BELLE to determine 
the quantum numbers of their hadronic final states. 

Combining the two classes of final states they arrive at: 

BR(£?~ — > ruD^n) = (1.81 ±0.20 ±0.20)% (50) 
BR(B d -> rPD^n) = (1.47 ±0.20 ±0.17)% (51) 

leaving room for a large D^irn component of ~ (1.3±0.4)%, whereas previous stud- 
ies have obtained 90% C.L. upper limits ranging from 0.35 to 1.3 %, see Eqs. (43,44). 

• BABAR's most recent analysis [6] yields rather consistent numbers: 



BR(B- - 


-> l~vDix) 


= (0.63±0.09 stat ±0.05 SJ/st )% 


(52) 


BR(B- - 


> rvD*n) 


= (0.89±0.08 stat ±0.06 SJ/st )% 


(53) 


BR(B d - 


-> VvDix) 


= (0.65±0.12 stat ±0.05 SJ/st )% 


(54) 


BR(B d - 


> rvD*ix) 


= (0.72±0.12 stat ±0.06 SJ/st )%, 


(55) 



which can be combined to 

BR(B~ -> rpD^n) = (1.52 ±0.12 ±0.10)% (56) 
BR(B d -> Z - PD W 7r) = (1.37 ± 0.17 ± 0.10)% (57) 

again leaving room for a significant D^nir component of about 1.3%. 

• Using DELPHI'S numbers stated in Eqs. (15, 45) and assuming that the "1" and "0" 
state decay 100 % into D*n and Dn, respectively, one arrives at 

BR(B- -> ruDn) = BR(B d -> I'DDn) ~ (0.9 ± 0.7)% (58) 
BR(B~ -> l-uD*n) = BR(B d — > l~9D*ir) ~ (1.9 ± 0.4)% . (59) 
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for a total of 



BR(£- -> ruD^n) = BR{B d -> Z _ PD W 7r) ~ 2.8% (60) 

One should note that the qualitative trend is the same as with BELLE's findings, 
Eqs.(46 - 49) - namely that D*tt final states dominate over Dn ones - yet the total 
D^it rate exceeds that reported by BELLE. 

The BT model predicts for Dn and D*n production: 

BR(B- -> Z-PDtt) = BR(B d -> Z _ P£>7r) = 0.51% (61) 
BR(5~ -> ruD*7i) = BR(B d -> ruD*7i) = 0.65% , (62) 

which is on the low side of BELLE's numbers and a forteriori for DELPHI'S findings, 
but consistent with BABAR's data. This is of course a rephrasing of the '1/2 vs. 
3/2' puzzle. 

In the BPS approximation [11] one has Ty\ = 0. Assuming that the sum rule of 
Eq.(18) saturates already with the n = state, one obtains Ty 2 = | leading to 

BR{B' -> ruDn) = BR(B d -> rpL>7r) = 0.39% (63) 
BR(B- -> rvD*7i) = BR(B d -> l~uD*Tr) = 0.50% ; (64) 

i.e., lower still. This might not be that surprising, since the BPS ansatz is at best 
an approximation rather than a systematic expansion. 

So far there is no experimental evidence for high mass hadronic states. One finds 
that about 6.4% and 18.3% of all D** states have masses between 2.6 and 3.3 GeV 
for the CDF and DELPHI data, respectively, which drop to 3.2 % and 7.8% for the 
mass range 2.8 to 3.3 GeV and 0.3 % and 3.1 % for 3.0 to 3.3 GeV. 

On the other hand, CDF seems to see more events below the D 3 ^ 2 peaks. Such low 
mass D^ir events could be due to higher mass states decaying into D^nn. CDF 
has not incorporated this scenario into their analysis, since previous measurements 
showed no evidence for such decays. 

One would conjecture that if the observed mass spectrum indeed differs significantly 
from theoretical expectations - in its center of gravity as well as its spread -, then 
the measured hadronic mass moments should not follow theoretical predictions - 
yet they do [2, 3, 4, 29, 33, 34, 36]. 

There is a more general problem. We have said before that the OPE treatment could 
a priori accommodate significant contributions with hadronic systems exhibiting a 
rather broad distribution in mass and centered below 2.5 or even 2.4 GeV. Yet 
closer scrutiny cast serious doubts on such a scenario. The OPE treatment involves 
applying quark-hadron duality. The latter can be expressed as saying that a rate 
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evaluated on the quark-gluon level can be equated with the sum of observable ex- 
clusive hadronic channels, at least after some averaging or 'smearing' over energy 
scales has been applied; the latter is sometimes referred to as semi-local duality. 

The question is: Which exclusive channels could be seen as dual to contributions 
with hadronic mass below 2.4 or even 2.5 GeV? (i) The theoretical estimates agree 
that the D 1 / 2 states that can populate this mass range possess too small production 
amplitudes to contribute significantly; their amplitudes actually would have to be 
enhanced greatly to overcome their kinematic suppression in semileptonic B decays, 
(ii) There are many higher orbital excitations of course, and taken together they 
might have a 'fighting' chance to yield a significant contribution to Tsl(B) - yet 
they all lie above 2.5 GeV in mass. It would not correspond to our usual picture 
that exclusive channels all above 2.5 GeV are dual to quark-level contributions 
computed to all lie below it. (iii) There is one intriguing possibility: We know of 
one basic failure of all quark models: they cannot explain the mass of the baryonic 
Roper resonance. For all quark models predict the mass of the first radial excitation 
higher than that for the first orbital excitation - an inequality clearly reversed for 
the Roper resonance. Would it be possible that the spectrum of charm resonances 
exhibits an analogous effect for mesons meaning that radial excitations can lie below 
the P wave states discussed before, and they make up the Dx component? This 
would be a most intriguing - yet also most exotic explanation. 

In summary: ALEPH, DELPHI and DO agree in finding a rate of about 0.8 - 1 % of Tb 
for the production of the two narrow D** states combined. This value is quite consistent 
with theoretical expectations for the D 3 ^ 2 rates. BELLE's data also fit naturally into 
this picture. The problem arises in the production of the broad D** states: The rate 
found by ALEPH and DELPHI suffice to saturate T SL (B), yet exceed the predictions 
for B — > WD 1 ! 2 by about an order of magnitude. BELLE's numbers on the other hand 
agree reasonably well with predictions, yet fall short of saturating Fsl(B). The mass 
distributions of the broad D** states, for which there is no clear experimental verdict, 
might pose a theoretical conundrum: if it is centered below 2.5 GeV, we have no natural 
candidates for these states. 

8 Comments on nonleptonic B decays 

It is usually argued - with very valid reasons - that the theoretical description is much 
murkier for nonleptonic than semileptonic B decays. We might encounter here one of the 
few exception to this general rule of thumb, and we have been alluding to this possibility 
already. One can analyze the inclusive transition 

B -> nX c (65) 

and study the hadronic system X c in the spirit of factorization; i.e., one analyzes its 
recoil mass spectrum, its quantum numbers and decay characteristics with the following 
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motivation: (i) The higher complexity of nonleptonic dynamics can be seen as an actually 
advantage here. For it provides additional production scenarios. More specifically it 
allows for significant production of D 1 / 2 states through the W emission diagram, which 
is not possible in semileptonic transitions, (ii) Without a neutrino in the final state it 
might be easier to determine the mass distributions and quantum numbers of X c . 

As already mentioned, for X c = D*, _D 3 / 2,1 / 2 the theoretical predictions from the BT 
model have been found in reasonable, though not compelling agreement with the data. 

One wants to extend such studies to the fully inclusive case and analyze also higher 
mass D** configuration. Intriguing first steps in this direction have already been under- 
taken by BABAR [22, 35]. The recoil mass spectra obtained by BABAR show clear D 
and D* peaks for charged and neutral B decays. The former exhibit also a clear peak 
around 2.4 - 2.5 GeV and maybe a signal also in the region above 2.6 GeV. In the Bd case 
there might be a signal in the 2.4 - 2.6 GeV region, but not much else. The peaks in the 
2.4 - 2.5 GeV region are natural candidates for showing D 3 ^ 2 production. The verdict on 
the domain beyond 2.5 GeV, which could be populated by the same configurations as in 
B — > l~vDfr oad , is tantalizing inconclusive. 

9 Conclusions and a call for action 

The Bd and B u inclusive semileptonic widths have been well measured. Most if not 
even all of it has been identified in B — > WD/D* + (0, l)ir. The theoretical description 
of B — > WX C rests on solid foundations. The potential discrepancies discussed in this 
note, which affect at most 20% of semileptonic B transitions, cannot lead to a significant 
increase in the uncertainty with which | V(cfo) | can be extracted from T(B — > WX C ). On 
the other hand they should not be 'brushed under the rug'. Theory does make non-trivial 
predictions of a rather sturdy nature. The OPE treatment is genuinely based on QCD, 
and while the BT description invokes a model, it implements QCD dynamics for heavy 
flavour hadrons to a remarkable degree. Their predictions therefore deserve to be taken 
seriously and not discarded at the first sign of phenomenological trouble. Preliminary 
lattice studies show no significant enhancement of '1/2' production. The numbers we 
have given for the theoretical expectations should be taken with quite a few grains of 
salt. Yet the predicted pattern that the abundance of '3/2' P wave resonances dominates 
over that for '1/2' states in semileptonic B decays is a robust one. Even a failure of such 
well-grounded predictions could teach us valuable lessons on non-perturbative dynamics 
and our control over them; it would certainly provide a valuable challenge to lattice QCD. 
Yet there is more: we know that D, D* and D 3 ! 2 production do not saturate T S l(B), and 
theory tells us that D 1 ! 2 cannot contribute significantly. What is then the nature of the 
missing hadronic configurations? 

On the experimental side the next important steps are: 

• In some of our discussion above we have modeled T SL (B) as the incoherent sum 
of B — > l~i)D / D* / D 3 / 2 / D x to infer the average mass of the configuration D x and 
its variance from the measured hadronic mass moments. The values of (M(D X )) 
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and y (M 2 (D X )) — (M(D X )) 2 serve as very useful diagnostics of the underlying 
dynamical situation. One finds that (M(Dx)) varies from 2.4 to 2.6 - or even 2.7 - 
GeV depending on whether one uses the branching ratios of Eq. (6) or of Eq. (5). 
The main reason for this relatively sizable shift is the variation in BK(B — > l~VD*). 
It would be most helpful to have this branching ratio determined more precisely. 

• Clarifying the size, mass distribution and quantum numbers of B — > lu[D/D*7r] b roa d 
and searching for B — > IvD / D* + 2ir with even higher sensitivity. 

• The data should be presented separately for B — > WD + 7r's and B — > IvD* + 7r's, 
since it provides more theoretical diagnostics. 

• More detailed analysis of B — > nX c , in particular in the high mass region for the 
X c and separately for charged and neutral B decays [22, 35]. 

These are challenging experimental tasks, yet highly rewarding ones as well: 

• They probe our theoretical control over QCD's nonperturbative dynamics in novel 
and sensitive ways. This is an area where different theoretical technologies - the 
OPE, quark models and lattice QCD - are making closer and closer contact. 

The lessons to be learnt will be very significant ones, no matter what the eventual 
experimental verdict will be: 

— A confirmation of the OPE expectations and even the more specific BT pre- 
dictions would reveal an even higher degree of theoretical control over nonper- 
turbative QCD dynamics than has been shown through T(B — > WX C ). 

— Otherwise we could infer that formally nonleading 1/mg corrections are highly 
significant numerically. Those corrections had to be highly enhanced to over- 
come the kinematic suppression in the production of D 1 ! 2 in semileptonic B 

decays. 

Such an insight would be surprising - yet important as well. In particular 
it would provide a highly nontrivial challenge to lattice QCD. Meeting this 
challenge successfully would provide lattice QCD with significantly enhanced 
validation. 

The call for further action is directed to theorists as well: 

• In the BT model one can compute the production rates for the higher orbital and 
radial excitations in semileptonic B decays. The individual rates seem to be rather 
small. It is conceivable that summing over a multitude of them yields a significant 
contribution. 

• The BT model predictions were obtained in the heavy quark limit. Corrections 
to this limit could be quite important as suggested in Ref. [24], and they could 
significantly change the relative weight of r[f 2 and r^\. Calculating or at least 
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constraining those corrections would be a most worthwhile undertaking - alas it 
requires some new ideas. A priori one can conceive of different ways of extending the 
BT description to include finite mass effects, yet they are unlikely to be equivalent. 
The foundations for a promising way have been laid in Ref.[15]. 

The authors of Refs. [37] find significant rates for the production of radial excitation 
that are enhanced further by 1/rriQ corrections. While we cannot quite follow their 
argumentation, 1/toq corrections due to potential exchange diagrams could indeed 
be sizable. This intriguing possibility needs and deserves intense scrutiny. 

There is also the exotic possibility that one finds here a mesonic analogue to the 
Roper resonance, namely that the radial (and even other) excitations of charm 
mesons are significantly lower in mass than predicted by quark models. This would 
presumably raise their production rates in semileptonic B decays to the sought-after 
level. 

• Lattice QCD studies of '1/2' and '3/2' production in semileptonic B decays has to 
be pursued with vigour. Such studies could turn out to be veritable 'gold mines' as 
far as validation is concerned. One can evaluate the spectrum of the higher radial 
and orbital excitations D**, for which some encouraging results have already been 
obtained [25] . Lattice calculations at finite values of m c should be performed, which 
would teach us about l/m c corrections. 

• The strong decays D** — > D/D* + mr should be estimated using heavy quark 
symmetry arguments augmented by quark model considerations. 

A final comment: The experimental analyses we advocate require considerable effort. 
We strongly belief such an effort is mandated by the insights to be gained, even if they 
are of a subtle nature: 

• They can provide us with important insights into the workings of nonperturbative 
dynamics. Lessons on the significance of I/tuq corrections and on systematic short 
comings of quark models would be of general value for the theoretical control we 
can establish over heavy flavour dynamics. 

• They probe the subtle concept of quark-hadron duality in novel ways. 

• The greatest practical gain might emerge for lattice QCD: If the latter can meet 
the challenge of such detailed data successfully, it would have gained a qualitatively 
new measure of validation. 

Acknowledgments: This work was supported by the NSF under grant number PHY- 
0355098 and by the EC contract HPRN-CT-2002-00311 (EURIDICE). 



19 



References 

[1] N.Uraltsev, hep-ph/0409125. 

[2] O. Buchmueller, H, Flaecher, hep-ph/0507253. 

[3] The Belle Collab., K. Abe et al, BELLE-CONF-0669. 

[4] The BABAR Collab.: B. Aubert et al, arXiv: 0707.2670 [hep-ex]. 

[5] HFAG, E. Barberio et al, hep-ex/0704.3575. 

[6] The BABAR Collab.: B. Aubert et al, arXiv: 0708.1738 [hep-ex]. 

[7] A. Czarnecki, M. Slusarczyk, F. Tkachov, hep-ph/0511004. 

[8] The BABAR Collab.: B. Aubert et al, Phys. Rev. D70 (2004) 091106. 

[9] I. Bigi et al, hep-ph/0512270. 

[10] I. Bigi, M.Shifman and N.G. Uraltsev, Ann. Rev. Nucl. Part. Set. 47 (1997) 591. 

[11] N.Uraltsev, Phys. Lett. B585 (2004) 253; ibid. B545 (2002) 337. 

[12] N. Uraltsev, Phys. Lett. B501 (2001) 86. 

[13] A. Le Yaouanc, L. Oliver, J.-C. Raynal, Phys. Rev. D67 (2003) 114009; Phys. Lett. 
B557 (2003) 207. 

[14] The BABAR Collab.: B. Aubert et al, hep-ex/0703027. 

[15] LI. Bigi, M. Shifman, N.G. Uraltsev and A. Vainshtein, Phys. Rev. D52 (1995) 196. 
[16] A. Oyanguren, talk given at EPS '05, Lisbon. 

[17] B. Bakamjian, L.H. Thomas, Phys. Rev. 92 (1953) 1300; A. Le Yaouanc, L. Oliver, 
O. Pene, J.-C. Raynal, Phys. Lett. B365 (1996) 319. 

[18] V. Morenas, A. Le Yaouanc, L. Oliver, O. Pene, J.-C. Raynal, Phys. Rev. D56 (1997) 
5668. 

[19] F. Jugeau, A. Le Yaouanc, L. Oliver, J.-C. Raynal, Phys. Rev. D72 (2005) 094010. 

[20] See also the 'covariant light-front approach' of H.-Y. Cheng, C.-K. Chua, C.-W. 
Hwang, Phys. Rev. D69 (2004) 074025. 

[21] The BELLE Collab., K. Abe et al, hep-ex/0412072. 

[22] The BABAR Collab., B. Aubert et al, Phys. Rev. D74 (2006) 111102. 



20 



[23] The BABAR Collab., B. Aubert et al, Phys. Rev. D74 (2006) 092004. 

[24] A. K. Leibovich, Z. Ligeti, I. W. Stewart and M. B. Wise, Phys. Rev. D57 (1998) 
308, [hep-ph/9705467]. 

[25] A.M. Green, J. Koponen, P. Pennanen and C. Michael, Eur. Phys. J. C28 (2003) 
79; J. Foley, A. O'Cais, M. Peardon, S. M. Ryan, Phys. Rev. D75 (2007) 094503. 

[26] S. Negishi, H. Matsufuru and T. Onogi, Prog. Theor. Phys. 117 (2007) 275. 

[27] D. Becirevic et al, Phys. Lett. B609 (2005) 298; B. Blossier, A. Le Yaouanc, 
V. Morenas and O. Pene, Phys. Lett. B632 (2006) 319; Erratum ibid 645 (2007) 
476. More details can be found in the thesis of B. Blossier (in French): it can be 
downloaded from the website http:/ /tel. archives-ouvertes.fr/. 

[28] The ALEPH Collab., Z. Phys. C 73 (1997) 601. 

[29] The DELPHI Collab., J. Abdallah et al, CERN-EP-PH/2005-015, hep-ex/0510024, 
accept, f. publ. in Eur. Phys. J. 

[30] The DELPHI Collab., P. Abreu et al, Phys. Lett. B475 (2000) 407. 

[31] The DO Collab., V.M. Abazov et al, Phys. Rev. Lett. 95 (2005) 171803. 

[32] The BELLE Collab., K. Abe et al, hep-ex/0507060. 

[33] The BABAR Collab.: B. Aubert et al, Phys. Rev. Lett. 93 (2004) 011803. 

[34] M. Battaglia et al, Phys. Lett. B556 (2003) 41. 

[35] A. Zghiche, LAPP-H-2007/01. 

[36] The CDF Collab., D. Acosta et al, Phys. Rev. D71 (2005) 051103; hep-ex/0502003. 

[37] D. Ebert, R. Faustov, V. Galkin, Phys. Rev. D61 (2000) 014016; ibidem D62 (2000) 
014032. 



21 



